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Abstract 

We investigate the possibility of probing high reheating temperature scenarios at the 
LHC, in supersymmetric models where the gravitino is the lightest supersymmetric par- 
ticle, and the stau is the next-to-lightest supersymmetric particle. In such scenarios, the 
big-bang nucleosynthesis and the gravitino abundance give a severe upper bound on the 
gluino mass. We find that, if the reheating temperature is ~ 10 8 GeV or higher, the 
scenarios can be tested at the LHC with an integrated luminosity of 0(1 fb _1 ) at y/s = 7 
TeV in most of the parameter space. 



1 Introduction 



Reheating temperature Tr is one of the most important physical quantities in the cos- 
mology, which represents the highest temperature of the universe at the beginning of the 
radiation-dominated epoch. From the viewpoint of cosmological observations, Tr is hardly 
restricted; the only constraint is Tr>5 MeV in order to realize a successful big-bang nu- 
cleosynthesis (BBN) [l]. 1 On the other hand, from the viewpoint of the particle physics, 
Tr has various implications. First of all, it is determined by the inflaton decay rate T^. If 
interactions between the inflaton and lighter particles are suppressed by the Planck scale 
M P = 2.4 x 10 18 GeV such as ~ m^/Mp, the reheating temperature becomes Tr ~ 
yT^Mp ~ 10 10 GeV(m^/10 13 GeV) 3 / 2 . In some inflation models, typical values of the infla- 
ton mass ~ lO^ 15 ) GeV provide T R ~ lfj( 4 ~ 13 ) GeV (cf. [4, 5]). Secondly, the baryon 
asymmetry of the universe must be generated between the inflationary era and the BBN 
epoch, and baryogenesis mechanisms give constraints on Tr. For instance, the thermal lep- 
togenesis scenario with heavy right-handed neutrinos [6] requires Tr > O(10 9 GeV) [7]. 
Thirdly, in supersymmetric models, the gravitinos are produced in the thermal plasma. 
Thus, the so-called gravitino problems impose severe upper bounds on the reheating tem- 
perature [8]. 

In this letter, we discuss the possibility of probing high reheating temperature scenarios 
at the LHC. 2 In particular, we consider the supersymmetric (SUSY) models with conserved 
R-parity, where the gravitino is the lightest supersymmetric particle, and the stau is the next- 
to-lightest supersymmetric particle. 3 We further assume that there is no entropy production 
after the inflation. 

In such a scenario, the BBN and the gravitino abundance give an upper bound on the 
gluino mass [12, 13, 14]. As the reheating temperature is larger, the mass bound becomes 
severer, and thus the production cross section of the SUSY events is enhanced. The essential 
feature of the events is the existence of heavy charged tracks of long-lived staus. Since the 
staus tend to have a high transverse momentum and a low velocity, the events are almost 
background free, and they can be easily discovered at the LHC. 

lr The reheating temperature may be probed with gravitational wave background in future experiments [2]. 
The BBN prediction on the 6 Li abundance, affected by long-lived charged particles, may also be used as a 
probe of the reheating temperature [3] . 

2 Possibilities of constraining T R at colliders in gravitino/axino LSP scenarios have been discussed in Ref. [9, 
10], though detailed collider studies have not been performed. 

3 LHC signatures of such a scenario have been studied in relation to cosmology from different points of 
view [11]. 
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In this framework, we will see that if the reheating temperature is ~ 10 GeV or higher, 
the scenarios can be tested even in the early stage of the LHC experiments, i.e. with an 
integrated luminosity of 0(1 fb _1 ) at y/s = 7 TeV in a wide parameter space of the SUSY 
Standard Model (SSM). Furthermore, we will find that a large part of the parameter region 
which is consistent with the thermal leptogenesis is accessible at O(10 fb _1 ) at \fs = 14 TeV. 

2 Cosmological Constraints 

In this section, we study the cosmological constraints on the present scenario. In Sec. 2.1, 
it is shown that the gravitino abundance gives upper bounds on the gaugino masses once a 
gravitino mass and a reheating temperature are provided. In Sec. 2.2, we discuss the BBN 
constraints on the long-lived stau and derive an upper bound on the gravitino mass for a 
given stau mass. Combining these two constraints, we obtain upper bounds on the gluino 
mass, depending on the stau mass and the reheating temperature. 

2.1 Gravitino abundance 

The present energy density of the gravitino must be smaller than the observed energy density 
of the dark matter. The gravitino energy density from the thermal scattering is given by [15, 
16, 17] 4 

x °' 14 (iiw) 2+0 

where m 3 / 2 is the gravitino mass, m~, m^, and m~ are the physical masses of the Bino, Wino, 
and gluino, respectively. We used the one-loop renormalization group equations to evolve 
the running masses up to the scale \x = Tr. Here, the numerical coefficients are evaluated 
at Tr = 10 s GeV. Logarithmic dependences on Tr is omitted for simplicity, but they are 
included in the following numerical calculation. The reheating temperature Tr is defined 
as Tr = (7T 2 <?* (TR)/90)~ 1 / 4 -\/r^Mp, where g* is the effective number of massless degrees of 
freedom. In this letter, mg and are taken to be free parameters. Thus, given the Bino 

4 Eq. (1) potentially includes an 0(1) uncertainty [16, 18]. In the present scenario, the gravitinos are also 
produced from the late-time stau decay. However, its contribution is negligible in the parameter region of 
our interest, and we neglect it in the following discussion. The gravitinos are also directly produced by the 
inflaton decay [5], but we do not include its contribution since it depends on the inflation models. 
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and Wino masses, we obtain upper bounds on the gluino mass from Eq. (1) by using the 
observed dark matter density flDAih 2 < 0.121 at 2a [19]. 

2.2 BBN constraints 

The BBN constraint is represented by the mass, the lifetime and the abundance of the stau. 
For a given mass and abundance of the stau, there is an upper bound on the stau lifetime, 
t~ ~ 487rMp77?|/ 2 m~ 5 . On the other hand, the abundance is determined by the stau mass. 
Thus, once we provide the mass of the stau, we obtain an upper bound on the gravitino mass. 

2.2.1 stau annihilation via electroweak processes 

The stau abundance is determined not only by the stau mass but also by its interaction with 
the other particles. In the SSM, the electroweak processes usually dominate the annihilation 
rate of the stau. The abundance below the freeze-out temperature is then estimated as [20, 
12] 5 



(2) 



Here and hereafter, we define Y~ = Y~ + + Y~- = (n~ +n~)/s, where n~ and s are the number 
density of the stau and the entropy density, respectively. Then, the BBN constraint leads 
to [8], 

< f 0.4 GeV - 10 GeV (lOOGeV < m~ < 450GeV) . , 

m 3 /2 „ | io GeV - 20 GeV (450GeV < m~ < lOOOGeV) ' [ } 

where the 6 Li constraint from the stau catalysis [22] gives the bound for m~ < 450GeV, while 
the deuterium constraint from the hadrodissociation due to the stau decay dominates for 
m~ > 450GeV [8]. 

Fig 1(a) shows the upper bound on the gluino mass for several reheating temperatures. 
Here, since we consider the stau NLSP scenarios, the region of m~ < m~ is excluded. Further- 
more, the Bino and Wino masses are set to be close to the stau mass, m~ = = 1.1 x m~, 
which provides a conservative bound, since this minimizes the Bino/ Wino contributions to 
the gravitino abundance [see Eq. (1)]. Note that if m^^) is too degenerate with m~, the 
stau abundance is enhanced by the B (W) decay, resulting in more stringent bounds on Tr 



5 Eq. (2) has an 0(10%) uncertainty which depends on details of the SSM parameters [20, 12]. We also 
neglected the correction from the Sommerfeld enhancement, which is also of the order 10% [21]. 
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(or m~) [20]. Although the gluino mass can be less than m~ and above m~, the region 
is very narrow, and thus it is discarded in Fig 1, for simplicity. 

We found that the gluino mass can be as large as about 2.4 TeV (2 TeV) for m~ < 
1 TeV (500 GeV) when the reheating temperature is Tr ~ 10 s GeV. For a higher reheating 
temperature, Tr>3(5) x 10 8 GeV, there are upper bounds on both the gluino mass and the 
stau mass, m~< 1 TeV (0.7 TeV) and m~ < 750 GeV (500 GeV). We also notice that there is 
no allowed region for Tr > 10 9 GeV. 

In the figure, the stau mass ranges up to 1 TeV. Although a very high reheating tempera- 
ture, e.g. Tr > 10 9 GeV, is not allowed in this mass range, there might exist such parameter 
regions if we consider a heavier stau. However, extrapolating the BBN bound in [8], and 
taking account of the constraint on the gravitino abundance from the stau decay, we do not 
find any parameter space which satisfies Tr>3 x 10 8 GeV for m~ > 1 TeV as long as the 
electroweak interactions dominate the stau annihilation. 



2.2.2 stau annihilation via enhanced Higgs channel 

The stau annihilation via the Higgs production channel may be enhanced by a large Higgsino 
mass parameter [i and tan/3, the ratio of the up- and down-type Higgs vacuum expectation 
values [23, 24]. Actually, the trilinear coupling of the lighter stau with the lightest Higgs is 

C = -A-tfnh, (4) 

where the coefficient is 

A~ h = --^(^tan/3 + ^ T )sin2# r + 5m ' 



2M W ^ M w 
-gzM z 



— - + sin 6y/ cos 9 T — sin 6\y sin L 



2 



(5) 



Here, g, gz, Mw, Mz, and 9\y are the Standard Model parameters. A T is the stau trilinear 
coupling with the down-type Higgs boson and T is the stau mixing angle. We assumed the 
decoupling limit of the heavy Higgs. When /xtan/3 is large, the first term in the right-hand 
side becomes enhanced. 

The stau relic abundance is obtained by solving the Boltzmann equation. It is approxi- 
mately given by 

^, 0xl „-. 5 (~!)(^), 
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where (av) is the thermally averaged annihilation cross section. When the electroweak process 
is dominant, we obtain the result in Eq. (2). On the other hand, as //tan/3 increases, the Higgs- 
and the top-pair production channel is enhanced. Then, the cross section becomes 

T T 

where Y t ~ 1 is the top Yukawa coupling, and and ft, are defined as 



, V 1 - r h an s f (1 - r t /2)Vl - r t 

fh= (l-r h /2? d{1 - rh) > ft= (l-r fe /4P ° {1 ~ ^ (8) 

with rh = m\/rn~ and rt = mf/m~. In the following, we use nit = 173 GeV and m/j = 
120 GeV. 

However, such a large trilinear coupling A~ h may suffer from a disastrous charge/color 
breaking (CCB) of the Standard Model [23]. Actually, the trilinear interaction A~ h can 
provide another vacuum apart from our vacuum in the stau-Higgs field space. As A~ h is 
enhanced, the CCB vacuum becomes deeper than our vacuum, the potential barrier between 
the two vacua is lowered, and hence our vacuum can decay into the CCB vacuum. 

In order to avoid the CCB, the lifetime of our vacuum is required to be longer than the 
age of the universe. The vacuum decay rate is evaluated by the bounce method [25] at the 
zero temperature, T = [23]. In addition, the vacuum can transit through thermal effects in 
the early universe. We analyzed the thermal transition by the method explored in Ref. [26]. 
The Higgs potential is evalulated at the one-loop level both for T = and T > 0. We also 
take account of the thermal potential, which comes from the top quark and the electroweak 
gauge bosons. As a result, we obtain an upper bound on A~ h as a function of m~. 6 Then, 
by using Eqs. (6) and (7), the minimum abundance of the stau is obtained, as is shown in 
Fig.2. In detail, the constraint comes from the vacuum decay at T = for m~>220GeV, 
while the finite temperature effect dominates for m~<220GeV. We approximated a bounce 
configuration by a straight line in the stau-Higgs field space. In fact, we checked at several 
model points that the approximation is accurate; the error is 0(1)% compared to the full 
two-dimensional result. 

We found that, although Y~ decreases more than one order of magnitude compared to 
the electroweak dominant case, since the stau abundance cannot be less than 3 x 1CT 15 for 
m~ > lOOGeV, the 6 Li overproduction from the stau catalysis still provides the severest 



6 The result is independent of the sign of fi, because if we flip it, the stau mixing angle 8 T changes the sign, 
too. 
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constraint. The Li bound depends only on the lifetime and the abundance of the stau, 
which can be read from Fig. 13 of Ref. [8]. By using the result in Fig. 2, we obtain 

™3/2 < 0.9GeV - 115GeV (lOOGeV < m~ < lOOOGeV). (9) 

Combined with the constraint from the gravitino abundance in Eq. (1), we obtain upper 
bounds on the gluino mass as shown in Fig. 1(b). Note that the gluino bound is similar to 
that of the electroweak case for m~ < 450GeV, because the constraint comes from 6 Li in 
both cases, while it is greatly ameliorated for m~ > 450GeV due to the absence of deuterium 
constraint. As a result, there is a region where a high reheating temperature, Tr > 10 9 GeV, 
is realized, which is favored by the thermal leptogenesis [23, 24]. 

3 Collider signatures 

In this section, we study the Tevatron constraints and LHC signatures of the present scenario. 
Since the stau is long-lived, the SUSY events predict charged tracks in the detectors. In the 
following analysis, we focus on the gaugino productions at the pp or pp collisions by assuming 
that the Higgsino and the scalar superparticles other than the stau are heavy enough not to 
contribute to the cross section. Then, the gauginos immediately decay into the stau in the 
end of the decay chains. As long as the squarks have masses of 0(1 — 10)TeV, the events 
are insensitive to the detail of the mass spectrum of the heavy superparticles. Furthermore, 
although the cosmological constraints in the last section depend on the size of A~ h , the 
production cross section of the stau is almost independent of it. 

For the numerical analysis, we use PYTHIA 6.4.22 [27] to simulate the kinematics of the 
events, and relied on PGS4 [28] for the detector simulation. The cross sections are estimated 
at the NLO level by Prospino 2.1 [29] for the gluino pair production, while the other (non- 
colored) channels are calculated at the LO by PYTHIA. 

3.1 Tevatron constraints 

The stable charged massive particles (CHAMPs), i.e. the stau in the present scenario, have 
been already constrained by the Tevatron experiments. Such a long-lived CHAMP looks like 
a heavy muon in the detector. In order for the particle to be distinguished from the muon, the 
experiments have searched for slow CHAMPs with a high transverse momentum. According 
to the CDF result [30], the events were selected by the following trigger: the highest px 
"muon" candidate has pr > 20 GeV, satisfying an isolation criteria Et (OA) /pt {muon) < 
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0.1, where J2^t(0A) is the sum of the transverse energy within a cone R = 0.4 around 
the candidate, excluding the energy which is deposited by the candidate itself. The CDF 
obtained the constraint that the production cross section of the non-colored CHAMPs with 
|r?| < 0.7, p T > 40 GeV and 0.4 < /3 < 0.9 is bounded to be less than 10 fb at the 95% C.L.. 

We show the CDF bound in Fig. 1. When the gluino is relatively light, the staus are 
produced via the gluino channel, while the neutralino/chargino productions become the main 
mode for a heavier gluino. In both channels, almost all the passed events are triggered by 
isolated staus. Here, we assumed = = 1.1 m~ (see Sec. 2.2). As and increase 
with m~ fixed, the CDF bound on the stau mass is weakened because the chargino/neutralino 
production cross section decreases, while the cosmological constraint on m~ (or Tr) becomes 
severer (see Eq. (1)). Although the bound is based on the CDF data with the integrated 
luminosity 1 fb" 1 , we checked that the sensitivity does not change greatly even for 10 fb -1 . 

3.2 LHC signatures 

Now let us discuss the LHC signatures. In the early stage, the LHC runs at the beam energy 
y/s = 7 TeV with an integrated luminosity up to about 1 fb -1 , and it is planned to be 
upgraded to ^fs = 14 TeV. In Fig. 1, we show the contours of the produced stau numbers at 
the LHC for 1 fb _1 at ^fs = 7 TeV and 10 fb" 1 at ^fs = 14 TeV. Similarly to the Tevatron, 
the vertical line comes dominantly from the pair production of the lightest charginos or 
that of the lightest chargino and the second lightest neutralino, while the horizontal line is 
obtained by the gluino production. As in the Tevatron case, if m~ and are increased 
from mg = = 1.1 m~, the vertical lines move to the left (i.e., less staus are produced), 
while the bound on m~ (or Tr) becomes severer. It should be mentioned that we imposed 
no cuts and triggers, which will be studied below. We find that much more staus will be 
produced even in the early stage of the LHC compared to the Tevatron. 

The realistic analysis includes the trigger criteria and the cuts to reduce the backgrounds. 
The particles produced in the decay chain can be energetic enough to trigger the events. 
Moreover, the stau itself may hit the muon trigger. In this letter, we adopt the following 
trigger menu, 

• at least one isolated electron has px > 20 GeV, 

• at least one isolated muon has px > 40 GeV, 

• at least one isolated tau has px > 100 GeV, 

• at least one isolated stau has px > 40 GeV within the bunch, 
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• at least two staus have pr > 40 GeV within the bunch. 

If any one of these conditions is satisfied, the event is read out. Here, the isolation condition 
relies on PGS4, but the muon and the stau are considered to be isolated if they satisfy the 
following two conditions: 1) the summed p^ in a R = 0.4 cone around the particle (excluding 
the particle itself) is less than 5 GeV, and 2) the ratio of £r in a 3 x 3 calorimeter array 
around the particle (including the particle's cell) to the pt of the particle is less than 0.1125. 
In the last two trigger items, we imposed the bunch condition, because the stau trigger works 
only when the stau reaches the muon trigger before the next bunch collides. Thus, the stau 
velocity is required to be f3 > 0.7 in the barrel region, \rj\ < 1.0, and (3 > 0.8 in the endcap 
region, 1.0 < \r/\ < 2.8, for the trigger to be applied [31]. Although we set relatively high 
thresholds for the triggers, the event number is almost insensitive to them. In addition to 
the above menu, one may take the jet triggers into account. We checked that this does not 
change the number of events for yfs = 7 TeV, while it can increase by an 0(1) factor for 
V« = 14 TeV. 

After selecting the events based on the above trigger menu, the stau number is counted 
with the following cut conditions for the stau: 

• pr > 20 GeV 

• 0.5 < p < 0.9 

• \r}\ < 2.5. 

Actually, the first two cuts significantly reject the muon background, because a high pt muon 
has a velocity of j3 ~ 1. The last one is imposed for the stau to be detected in the muon 
detector [31, 32]. 

The number of events may be further reduced by taking the efficiency into account. In 
the current setup, a large part of the events passes the muon (stau) triggers. In fact, the 
most events are triggered by the stau when the chargino/neutralino is the dominant channel. 
In the gluino production case, roughly a half of the total passed events is accepted by the 
stau trigger, and the others are by the leptons. Since the trigger efficiency of the muon-like 
event is pretty good [31], the reduction is expected to be small. On the other hand, the 
reconstruction efficiency of the stau depends on the LHC analysis. According to the ATLAS 
CSC studies [31], the efficiency varies from 0.1 to 0.9 for (5 ~ 0.5 — 0.9. The analysis can 
be improved by the method studied in [32], which provides the efficiency more than 90% for 
f3 > 0.5. In the following result we assume the efficiency to be 100% for simplicity. 
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In Fig. 3, we show the number of the staus for the LHC beam energy y/s = 7 TeV and 
an integrated luminosity 1 fb _1 . The result after taking account of the detector simulation 
is compared with that without any cut conditions. We find that the number decreases down 
to ~ 50% and 20 — 30% for Fig. 3(a) and (b), respectively. More precisely, the efficiency 
depends on details of the superparticle mass spectrum. 

The long-lived CHAMP events are almost free from the background, and therefore they 
can be discovered with a small number of the events. For instance, we see that the LHC will 
produce more than 5 (10) staus for m~<460GeV (410GeV) when the chargino/neutralino 
production is dominant (Fig. 3(a)), and the same for m~<990GeV (910GeV) in the case 
of the gluino production (Fig. 3(b)). In a smaller mass region, the stau number depends 
on both the stau and gluino masses. Especially, we checked that the whole region with 
Tr>3 x 10 8 GeV in Fig. 1(a) predicts more than 10 staus after applying the cut conditions. 
Furthermore, when the LHC beam energy is upgraded to \/s = 14 TeV and the integrated 
luminosity reaches 10 fb -1 , the sensitivity becomes greatly improved. Actually, we checked 
that more than 10 staus can be detected even for m~ = ITeV and m~ = 2.5TeV. 

In the above analysis, we assumed that the scalar superparticles, particularly the squarks, 
are very heavy. As the squarks become lighter, productions of the squarks come into the chan- 
nels. In fact, when the squark has a mass comparable to the gluino mass, the squark channels 
become dominant, and the SUSY cross section is enhanced by one order of magnitude com- 
pared to that of the gluino. Thus, it is considered that the result in this letter is conservative, 
and we may observe more SUSY events even in the early stage of the LHC. 

4 Conclusion and Discussions 

In this letter, we studied the cosmology and the collider features of high reheating temper- 
ature scenarios with a long-lived stau within the SSM. We first discussed the cosmological 
constraint from the gravitino abundance and the BBN. In most of the parameter space, the 
stau annihilation is dominated by the electroweak interactions, and we obtained the upper 
bounds on the gluino mass for a high reheating temperature. When the reheating tempera- 
ture is Tr > 10 8 GeV, the gluino is limited to be lighter than about 2.4 TeV for m~ < ITeV. 
A higher reheating temperature provides a stronger constraint, m~< \ TeV (0.7 TeV) and 
m~< 750 GeV (500 GeV) for T R >3(5) x 10 8 GeV In addition, there is no allowed region for 
Tr > 10 9 GeV. We then investigated the collider signatures. We found that more than 10 
staus will be detected for Tr>3 x 10 8 GeV in the early LHC for an integrated luminosity 
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1 fb 1 with y/s = 7 TeV. Upgrading the LHC to <fs = 14 TeV, we expect > 10 staus in the 
whole mass region of Tr> 10 8 GeV with m~ < ITeV for an integrated luminosity 10 fb -1 . 

When the Higgs coupling to the stau is very large, the cosmological constraint becomes 
ameliorated, especially for a heavy stau. Interestingly, this opens a window for a high re- 
heating temperature, TR>10 9 GeV for m~>400GeV, which saves the thermal leptogene- 
sis [23, 24]. We found that this region predicts more than 5 staus in the early LHC if the stau 
mass is m~ < 700GeV, and the stau number becomes much more than 10 for m~ < ITeV till 
the end of the LHC. Thus, the thermal leptogenesis in the present framework can be tested 
by the LHC in a wide parameter region of the SSM. 

If the long-lived staus are indeed produced at the LHC, the masses of the stau and other 
superparticles can be measured with a good accuracy [33] . It then becomes possible to obtain 
a more precise (severer) bound on the reheating temperature. Furthermore, the stau lifetime 
(and hence the gravitino mass) may also be measured [34]. The stau lifetime is directly 
subject to the BBN constraint, while the gravitino mass gives a further information on the 
gravitino abundance. In particular, if there is no allowed region in the parameter space of 
the stau lifetime and the stau abundance, it indicates a very low reheating temperature [3] 
or an entropy production between stau decoupling and the BBN [35]. 

Although we studied the electroweak and the Higgs production channels of the stau 
annihilation, it is difficult to distinguish them merely by counting the stau number at the 
LHC. This is because the event kinematics is almost insensitive to the stau-Higgs interaction. 
Noting that the lighter stau is likely to mix with both the left- and right-handed staus 
maximally in order to realize a large A~ h , we may measure the mixing angle by determining 
the chirality structure of the stau, for example, by measuring a polarization of the tau which 
is associated with the stau production [36]. 

In this letter, we discussed the stau annihilation via the electroweak interactions and the 
light Higgs channel. Ref. [24] also studied a possible annihilation channel of the stau through 
the heavy Higgs exchange. If the masses are tuned well such that the stau annihilates through 
the heavy Higgs pole, the stau abundance may decrease considerably. This scenario may be 
interesting from the viewpoint of the LHC signatures, since the heavy Higgs has a strong 
interaction with the stau and its mass is correlated with the stau mass. 

When the stau coupling to the Higgs is strong, the SUSY contribution to the muon 
anomalous magnetic moment may become large. This is because large \x and tan/3 enhance 
the mixing of the left- and right-handed sleptons, namely the smuon-neutralino contribution 
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to the muon g — 2 [37]. Since the experimental result [38] shows a 3-4cr discrepancy from the 
SM prediction [39], a large A~ h might be an interesting possibility. 
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Figure 1: The solid (blue) lines are upper bounds on the gluino mass for various reheating 
temperatures Tr with m~ = = 1.1 m~. The stau annihilation is dominated by (a) 
the electroweak interactions and (b) the Higgs production channel. On the other hand, the 
dotted (black) and the two-dotted dashed (green) lines are contours of the produced stau 
numbers at the LHC with the integrated luminosity 1 fb _1 at y/s = 7 TeV and 10 fb _1 
at y/s = 14 TeV, respectively. Here, the SUSY events only include the productions of the 
gluino, chargino and/or neutralino (and stau), and no cuts and triggers are imposed. The 
dotted-dashed (red) line is obtained from the CDF bound. Under the dashed (purple) line, 
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Figure 2: The stau abundance Y~ for varying the stau mass m~. The electroweak interaction 
gives the stau annihilation on the dashed line, while the Higgs (and top) production channel 
is dominant on the solid line. Here, the stau coupling to the Higgs boson, A~ h , is taken to 
be maximal for the latter case. The top and Higgs masses are set to be m t = 173 GeV and 
m h = 120 GeV. 
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Figure 3: The solid line shows the number of the generated staus after applying the triggers 
and the cuts, which is compared to the stau number without any conditions displayed by the 
dashed line. The gluino mass is set to be 2TeV for (a), and the stau mass is 600GeV for (b). 
The beam energy is chosen to be yfs = 7 TeV. 
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